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ABSTRACT 

We consider the limiting case for orbital stability of the companions to HR 8799. This case is 
only consistent with ages for the system of ~ 100 Myr, not with the 1 Gyr age proposed from 
astroseismology. The discrepancy probably arises because the inclination of the star is smaller than 
assumed in analyzing the astroseismology data. Given this young age, the best estimates of the 
companion masses place them by a small margin on the planet side of the division between planets 
and brown dwarfs. 

Subject headings: circumstellar matter — infrared: stars — planetary systems — stars: HR 8799 



1. INTRODUCTION 

HR 8799 is a late A/early F star located at 39.4 pc, 
with T off ~ 7350K, M* = 1.5 M and a A Bootis-like pat- 
tern of surface metallicity. The star has received a very 
high level of attention since the imaging of three massive 
substellar objects in large radius orbits around it (Marois 
et al. 2008). HR 8799 and Fomalhaut (Kalas et al. 2008) 
are the only two firm examples of systems with separa- 
tions > 50 AU between the star and planetary-mass com- 
panions. Both stars also have luminous planetary debris 
disks that can provide further insight to the outer struc- 
tures of their planetary systems (e.g., Stapelfeldt et al. 
2004; Su et al. 2009). Given the bias in current planet- 
discovery techniques toward close companions, the exis- 
tence and properties of these widely separated examples 
open a new perspective on the diversity of planetary sys- 
tems. 

The nature of the companions of HR 8799 depends 
critically on the age of the system, since their cooling 
rates (Baraffe et al. 2003) enter strongly into the relation 
of their luminosities to their masses. The star is near the 
zero-age main sequence, where traditional methods for 
estimating its age are not accurate. Marois et al. (2008) 
used a variety of arguments 4 to estimate an age of a 30 — 
160 Myr and on this basis assign masses of 5 — 11, 7— 13, 
and 7—13 Mj up respectively to HR 8799 b, c, and d. 
Janson et al. (2010) used VLT/NACO to obtain spatially 
resolved spectroscopy of HR 8799 c from 3.88 to 4.10 /jm. 
Although there are some discrepancies in the slope of the 
spectrum beyond 4 /iin, it can generally be fitted by a 
COND model from Baraffe et al. (2003) with T pl = 1100 
K, log g = 4.0 and R p ; = 1.3 R,/ up . These characteristics 
are consistent with the mass estimate for planet c of 10 
M Jup at 60 Myr by Marois et al. (2008). The minimum 
mass for deuterium burning, ~ 13 M,/ up , is often taken as 
the dividing line between planets and brown dwarfs, in 
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which case these mass estimates indicate that HR 8799 
hosts a fascinating planetary system. 

However, Moya et al. (2010a) argue that the stellar 
age is still unconstrained, pointing out that some of the 
age estimators are compromised by the A Bootis char- 
acter of the star or are statistical in nature. They used 
astroseismology to show that an older age of ~ 1 Gyr is 
possible, if the inclination of the equatorial plane of the 
star is i» > 36° . Because the masses of the companions 
are derived from the comparison of the observed luminos- 
ity to theoretical evolutionary tracks of giant planets and 
brown dwarfs, the age of the star is critical in determin- 
ing these masses. At this older age, the companion lumi- 
nosities would have faded significantly and their assigned 
masses would be well into the brown dwarf range. Thus, 
the nature of the HR 8799 system would change signifi- 
cantly. This would affect the formation scenarios under 
discussion, and the ability to use these three planets (of 
similar age and metallicity) to calibrate evolutionary and 
atmospheric models of giant planets. 

In this letter, we assess the possibility of more massive 
companions by considering the limiting case in which the 
masses of the planets are set to the largest values that 
dynamical stability can allow (Section [2]). Given the 
observed planetary luminosities and using evolutionary 
models, we estimate the age corresponding to these lim- 
iting masses. We then evaluate the stability of the plan- 
ets and the dust-producing planetesimals during that 
timescale to assess whether an older age could still lead 
to a planetary system (planets + planetesimals) consis- 
tent with the observations (Section [3]). We discuss the 
results in Section 21 

2. A STABLE CONFIGURATION WITH MASSIVE BODIES - 
A LIMITING CASE 

We now consider whether the configuration of the com- 
panions can help resolve the uncertainty over the age 
of HR 8799. The companions have been detected in 
HST/NICMOS archival data taken 10 years prior to the 
discovery image (Lafreniere et al. 2009), however, their 
long orbital periods make the determination of the or- 
bital parameters still uncertain. While there are several 
configurations that would fit (1) the astrometry, (2) the 
observed luminosities, and (3) the requirement that the 
system is stable for at least 30-160 Myr, dynamical anal- 
ysis by various authors converge on a solution in which 
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TABLE 1 

Planetary and stellar parameters for the dynamical 

SIMULATIONS (FlT D6A) 



Mplanet 

(M ) 


a 

(AU) 


e 


i 

(rad) 


n 

(rad) 


w 
(rad) 


M 
(rad) 


1.27-10" 2 (b) 


67.91 


0.002 








7T 


7T 


1.80-10" 2 (c) 


37.97 


0.005 














1.80-10" 2 (d) 


23.52 


0.083 








7T 






Note. — Planetary parameters used in the dynamical sim- 
ulation D6a (taken from Fabrycky & Murray-Clay 2010). The 
planet masses correspond to 13.30, 18.84 and 18.84 Mj„ p , for 
planets b, c and d, respectively, a and e are the semi-major 
axis and eccentricity of the planet; the orbits are face on (I p i 
= 0) and co-planar (i = 0); ui is the argument of periastron, 
Q is the longitude of the ascending node, and M is the mean 
anomaly. The mass of the central stars is 1.5 Mq. 

the companions are on nearly circular orbits, locked into 
a double 4:2:1 mean motion resonance (ld:2c:4b). This 
possibility was first suggested by Fabryck & Murray- 
Clay (2010), and latter confirmed by Gozdziewski & Mi- 
gaszewski (2009) and Reidemeister et al. (2009). (For 
a summary of the possible orbital solutions see section 
2.4.2 in Moro-Martin et al. 2010). The companions avoid 
close encounters in the resonant configuration, ensuring 
dynamical stability. 

In this paper, we focus on this configuration because it 
may be the only solution that can guarantee the dynam- 
ical stability of the system on a timescale comparable to 
the stellar age, and in the case of massive planets. The 
analysis in Fabrycky & Murray-Clay (2010) concluded 
that, under this configuration, the planetary masses can 
be increased up to 1.9 times their nominal values (see 
their Figure 11). Therefore, we consider the following 
planetary masses with orbital elements listed in Table [TJ 
1.27-10" 2 Mq (= 13.30 M Jup ) for planet b, and 1.80-10~ 2 
M© (= 18.85 Mj up ) for planets c and d; this would put 
planets c and d well into the brown dwarf regime, while 
planet b would be near the transition. Following the la- 
beling convection in Moro-Martin et al. (2010), we will 
refer to this configuration as fit D6. 

The conclusion that the values adopted here are the 
maximum permitted masses is based on the analysis by 
Beauge et al. (2003): for a given pair of eccentricities 
(ei, e2), they estimated the maximum mass of the inner- 
most planet (Mi) for which stable apsidal corotation at 
the 2:1 resonance can exist. For the configuration consid- 
ered in this paper (see Table [JJ), where the planet pair b-c 
has ei = e c = 0.005 and e-i = e\> = 0.002, Beauge et al. 
(2003) estimated that the maximum mass of Mi — M c is 
< 0.012 M„ (see their Figure 8); similarly, the planet pair 
c-d, with e\ = ed = 0.083 and e% = e& = 0.005, would 
correspond to M\ = Md having a maximum mass of ~ 
0.013 M*. These values are similar to those adopted in 
this paper for planets c and d (1.80 -10~ 2 M = 0.012 M*, 
for M*= 1.5 M ). Configurations with planetary masses 
more massive than these cannot have a stable corota- 
tion orbit in the 2:1 resonance. In other words, the fit 
D6 adopted in this study is a limiting case: the upper 
bound for the dynamically allowed planetary masses. 

Using different evolutionary models, and based on the 
luminosities reported by Marois et al. (2008), log(-^) 

= -5.1 for planet b and log(-j^) = -4.7 for planets c and 



d, we can estimate the age of the system that would 
correspond to the observed luminosities and the masses 
adopted in this paper: using the models of Chabrier et 
al. (2000), the age would be 100-500 Myr (it cannot 
be constrained any further because their tables do not 
include results for intermediate ages); from Baraffe et 
al. (2003 - their Figure 1), we obtain an age of 250-350 
Myr; while from Burrows et al. (1997 - their Figure 7) 
we get that if planet b were to have 1.27T0 -2 M Q , its 
observed luminosity would correspond to a stellar age of 
~ 260 Myr. The temperature of 1100 K for HR 8799 c 
(Janson et al. 2010) can also be used to constrain the 
age: from Chabrier et al. (2000 - their Figure 2), this 
temperature and adopted mass would correspond to an 
age of < 225 Myr, while from Baraffe et al. (2003 - their 
Figure 1) the age would be in the range 250-280 Myr. 
We conclude that if the planets were to have the limiting 
masses adopted in this paper (1.9 times their nominal 
values), their current luminosities would imply a system 
age of approximately 250-350 Myr. 

3. STABILITY OF COMPANION AND DUST-PRODUCING 
PLANETESIMAL ORBITS 

We now evaluate the stability of the planets and 
the dust-producing planetesimals on the timescales es- 
timated above. Following Fabrycky & Murray-Clay 
(2010), and to assess how sensitive the results are to small 
changes in initial conditions, we consider several plane- 
tary configurations around the one described in Section 
[2] and Table [TJ the differences are 10~ 4 AU for the semi- 
major axes, 10 -4 for the eccentricities, and 0.01° for the 
inclinations, ascending nodes, longitudes of periastron, 
and mean anomalies. In total, we model 14 different sets 
of initial conditions for the orbits. In addition to the 
three companions and the central mass of 1.5 Mq, the 
simulations include 500 test particles. Seven of the simu- 
lations consider particles in the inner disk, uniformly dis- 
tributed between 2 AU and the semi-major axis of the 
inner- most planet, while the remaining seven consider 
particles in the outer disk, spaced uniformly between the 
semi-major axis of the outermost planet and 300 AU. We 
assume that the planets and the dust-producing planetes- 
imals formed out of a thin disk and are co-planar. The 
particles are on initially circular orbits, with angular ele- 
ments chosen randomly between and 2tt. Particles are 
removed if they approach the star closer than 2 AU, or 
approach a planet closer than its Hill radius. The orbits 
were integrated using the multiple time step symplectic 
method skeel-SyMBA (Duncan, Levison & Lee 1998). 

3.1. Planets 

Figures [TJ and [2] show the evolution of the planets' 
semi-major axes, eccentricities, and resonant angles for 
the 2:1 resonance between planets c and d (4>d), and the 
2:1 resonance between planets b and c (0 C ), for the two 
longest-lived simulations, one for the inner disk and one 
for the outer disk, with dynamical lifetimes of 137.5 Myr 
and 147.5 Myr, respectively. While the companions are 
locked into the double resonance, the system is stable. 
Table [2] lists the resonant angles (mean and standard de- 
viation) and lifetimes of the 14 simulations ran, the lat- 
ter ranging from of 9 to 147 Myr. Using the MERCURY 
integrator, Fabrycky & Murray-Clay (2010) found that 
the time to instability of their six simulations (with the 
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Fit D6a (inner disk) 



80 




Time (Myr) initial semi-major axis (AU) 



Fig. 1. — Left: Long-term evolution of the planets' orbits un- 
der simulation D6a, during the time in which the system is stable. 
The orbital elements are listed in Table[l] From top to bottom, the 
panels correspond to semi- major axis (a), eccentricity (e) and res- 
onant angles for the 2:1 resonance between planets c and d (<j>d), 
and the 2:1 resonance between planets b and c (</>c). The mean 
and standard deviation of the resonant angles are listed in Table [2] 
The colors correspond to planet b (black), c (orange) and d (blue). 
Right: Results from the dynamical simulation of 500 test parti- 
cles in the inner planetesimal disk. Top: test particle's lifetimes. 
Middle: allowed parameter space for the planetesimals' orbital el- 
ements, where the shaded areas indicate regions where test par- 
ticle's orbits are unstable due to planet-crossing (striped area) or 
overlapping first order mean motion resonances (dotted area); the 
squared symbols show the maximum eccentricity attained by test 
particles on initially circular orbits. Bottom: number of surviving 
test particles at the onset of dynamical instability. 

Fit D6i (outer disk) 



80 




Time (Myr) initial semi— major axie (AU) 

Fig. 2. — Same as Figure [T] but for a simulation of 500 test 
particles in the outer planetesimal disk. In this case, the high 
baseline value of the test particles' eccentricities and slight rising 
trend is a numerical artifact from using stellar-centric osculating 
elements. 



same set of initial conditions) range from 25.9 Myr to 
102.6 Myr. We conclude that if the planets were to have 
the limiting masses adopted in this paper (1.9 times their 
nominal values), the lifetime of the system would be < 
150 Myr. 



3.2. Dust-producing planetesimals 

Following Moro-Martm et al. (2007, 2010), we have 
estimated the potential location of the dust-producing 
planetesimal belts under the configurations considered 
in this paper. The possible niches are the regions where 
the orbits are stable and the maximum eccentricities ex- 
cited by the planetary perturbations remain low (e max < 
0.3), ensuring the survival of planetesimals against quick 
erosion from mutual collisions. Figures Q] and [5] show the 
results for the two longest-lived simulations already dis- 
cussed (one for the inner disk, and one for the outer disk). 
We find that the outer edge of the inner (warm) planetes- 
imal disk is Rout" 1 ~ 9 AU, an d that the inner edge of 
the outer (cold) planetesimal disk is Rf° ld ~ 150 AU. 
We have adjusted the disk model presented in Su et al. 
(2009) to this new geometry and find that a satisfactory 
fit to the spectral energy distribution can be achieved 
readily. (For comparison, the configuration adopting the 
nominal masses of 7, 10 and 10 Mj up for planets b, c 
and d, respectively - labeled Dl in Moro-Martm et al. 
2010 - results in R™ rm - 12 AU and Rf° ld - 110 AU). 
We conclude that the presence of planets with masses 
1.9 times their nominal values would still be consistent 
with the SED and imagining of the debris disk around 
HR 8799. 

4. DISCUSSION AND CONCLUSION 

The configuration considered in this paper (D6) is a 
limiting case, in the sense that the adopted companion 
masses are the maximum values for which the stabiliza- 
tion mechanism, the double 4:2:1 resonance, can apply. 
Based on stability considerations, it is therefore unlikely 
that the companions are more massive than considered 
here. 

In Moro-Martm et al. (2010) we found that configu- 
rations in which the planets adopt their nominal masses 
(7, 10 and 10 Mj up for planets b, c and d, respectively - 
labeled Dl, D4 and D5) are consistent with planet and 
dust observations because: (1) the planets masses corre- 
spond to an age estimate of ~ 60 Myr and the system 
is found to be stable for > 160 Myr; and (2) the orbits 
of the dust-producing planetesimals fulfill two require- 
ments: (2a) they are in agreement with the dust location 
as derived from debris disk observations; and (2b) for 
the timescale during which the planet orbits are stable, 
the planetesimal orbits are also stable and their eccen- 
tricities remain low (implying survival against erosion). 
In contrast, the configuration considered in this paper, 
where the planets are 1.9 times more massive, is not fully 
consistent. Conditions (2a) and (2b) above are fulfilled: 
for the timescale during which the planets are on stable 
orbits, there are planetesimals in the regions where the 
dust is inferred to be located, and the planetesimal orbits 
are stable and their eccentricities remain low. However, 
condition (1) is violated: while the substellar object evo- 
lutionary models predict an age of the system around 
250-350 Myr (Baraffe et al. 2003) or - 260 Myr (Bur- 
rows et al. 1997) from luminosity constraints, or < 225 
Myr (Chabrier et al. 2000) and 250-280 Myr (Baraffe 
et al. 2003) from temperature constrains, the dynamical 
lifetime of the system is significantly smaller, < 150 Myr. 
Based on dynamical considerations we therefore conclude 
that the large companion masses considered in this paper 
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TABLE 2 

Results from dynamical simulations 



Model 

Dyn. Lifetime (Myr) 
mean(</> c ) (°) 

<K*c) (°) 

mean(^) (°) 



D6a 


D6b 


D6c 


D6d 


D6e 


D6f 


D6g 


D6h 


D6i 


D6j 


D6k 


D61 


D6m 


D6n 


137.5 


41.2 


110.9 


59.3 


29.9 


86.5 


32.0 


117.3 


147.5 


43.3 


19.2 


9.1 


17.6 


31.8 


-0.04 


-0.84 


0.14 


-0.50 


6.52 


-0.06 


-4.82 


0.86 


0.19 


-1.64 


1.03 


0.29 


0.67 


0.56 


44.4 


45.7 


44.6 


43.9 


48.7 


44.1 


48.9 


44.3 


43.4 


46.4 


103.2 


103.5 


103.1 


102.5 


-0.08 


-0.17 


-0.20 


-0.08 


1.90 


0.03 


-0.96 


0.24 


0.01 


-0.13 


-5.67 


3.47 


0.80 


-0.56 


8.9 


9.9 


9.2 


9.4 


10.3 


9.0 


9.9 


9.2 


8.6 


10.1 


103.0 


104.2 


101.0 


102.0 



Note. — Results for the dynamical simulations around the configuration described in Table [T] Dyn. Lifetime is the time for 
instability in Myr. Mean(cf> c ) and mean(0 ( j) are the mean of the resonant angles for the 2:1 resonance between planets b and 
c (4> c ), and the 2:1 between planets c and d ((f>d), during the time in which the system is stable. rf> c = 2A;, — A c — ui c and 
<f>d, = 2A C — A^ — [id, where A is the mean longitude, A = M + f2 + U3, and Co is the longitude of pericenter, uj = Q + u>. cr(<p c ) and 
cr(4>d) are the standard deviations of the resonant angles. The test particles in models D6a-e, k and m are initially located in the 
inner disk, while the test particles in models D6f-j, 1 and n are located in the outer disk. 



are unlikely and that the dynamical state of the system 
favors ages of ~ 100 Myr and companion masses slightly 
below the brown dwarf regime. 

This conclusion is supported by the strong emission of 
the debris disk at 24 /*m, behavior that is very uncommon 
around stars older than ~ 500 Myr (Caspar et al. 2009). 
In addition, the average temperature of 7350 K and the 
luminosity of 5 L from Sadakane (2006) and Moya et 
al. (2010b) places the star at or below the zero age main 
sequence and at a very young age compared with other 
A Boo stars (Paunzen et al. 2002). 

However, these results are not consistent with the pos- 
sibility that the system is ~ 1 Gyr old. The astroseismol- 
ogy modeling that suggests this age (Moya et al. 2010a) 
depends critically on the rotation speed of the star and 
hence on the assigned inclination. It is only valid if the 
inclination of the equatorial plane of the star is /* >36°. 
The case for /* = 36° results in two sets of possible so- 
lutions, a stellar mass of 1.32-1.33 M Q and an age of 
1123-1625 Myr, or a stellar mass of 1.44-1.45 M and 
and age of 26-430 Myr (with 16.7% of the models lying 
in the generally adopted age range of 30-160 Myr). For 
/* = 50°, Moya et al. (2010a) obtain a stellar mass of 
1.32 M and and age of 1126-1486 Myr. 

However, Lafreniere et al. (2009) suggested that the 
inclination of the orbital plane of the planets is I p \ = 13- 
23° from astrometry observations with a 10 year baseline 
(using NICMOS archival data). Su et al. (2009) reported 



from the degree of symmetry of the spatially resolved 70 
/im disk that the inclination of the disk is consistent with 
being face-on and unlikely to be /disk > 25°. We have re- 
evaluated these observations and place a 3a upper limit 
to id isk of 40°. Upcoming JCMT/SCUBA-2, Herschel 
and HST imaging of HR 8799 will be able to further 
constrain the inclination of the debris disk. If we as- 
sume that the equatorial plane of the star coincides with 
the orbital plane of the planets and the dust-producing 
planetesimals (i.e. i* = I p \ = /disk), these observational 
constraints indicate that indeed the inclination of the 
star is likely inside the 18 ° < /* < 36° range where the 
astroseismology models cannot be applied, or that it is 
close to 36° where they permit solutions consistent with 
the ~ 100 Myr age derived from other considerations. 

We conclude that HR 8799 is young, with an age of 
~ 100 Myr. Its companions are therefore likely all to 
be very massive planets, just below the mass limit for 
deuterium burning. 
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